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ABSTRACT 

The  hydrodynamic  force  and  yaw  moment  acting  on  a  slender 
axi symmetric  body  are  found  for  the  case  in  which  the  body  is 
moving  along  its  axis  witn  constant  forward  velocity  parallel 
to  a  vertical  wall.  The  fluid  is  assumed  to  be  inviscid  and 
incompressible.  7/.  ..  ■  Cr>’  •  c  f'/rs  ^  N 

First,  -the  situation  of  a  body  near  a  flat  wall  <£a-  addres-sed 
using  the  method  of  images.  The  body  and  its  image  are  modeled 
using  a  continuous  line  distribution  of  sources  and  doublets; 
the  dipole  distribution  is  used  to  account  for  velocities 
induced  on  each  body  by  its  image. 

Next,  the  case  of  a  body  running  parallel  to  the  mean  position 
of  a  wall  varying  sinusoidally  in  the  longitudinal  direction  is 
analyzed.  The_j._nteract ion  between  the  sinusoidal  wall  and  body 
is  modeled  using  a  "large"  axisymnetric  body  in  proximity  of  a 
smaller  oody.'i  The  quasi-static  case  is  investigated,  that  is, 
the  body  and  tfall  are  held  fixed  in  a  uniform  stream. 

The  transverse  force  and  moment  on  a  body  near  each  type  of 
wall  are  determined  us^ng  two  different  methods:  Lagally’s 
theorem  and  ^6egmenteda  t heory . y ( "Segment ed "  refers  to  dividing 
the  jody  into  vertical  segments  and  calculating  the  force  on 
each  segment  using  either  a  two  or  three-dimensional  flow 
analysis,  whichever  is  appropriate.  In  the  two-dimensional 
case,  the  segments  correspond  to  "strips"  used  in  strip 
theory .)  Lagally' s  theorem  is  used  to  find  the  axial  force  on  a 
ooay  near  a  sinusoidal  wall. 

Computer  programs  and  calculated  results  are  presented  for  an 
unappended  modern  submarine  hull  form  in  the  vicinity  of  both 
types  of  walls. \  Force  and  moment  are  found  to  increase  rapidly 
as  the  distance' between  wall  and  body  decreases.  In  the  flat 
wall  case,  Lagally  and  segmented  theory  calculations  correlate 
well  with  model  iest  results.  For  the  sinusoidal  wall  problem, 
results  are  plotted  indicating  how  force  and  moment  vary  with 
wall  amplitude  end  longitudinal  location  of  the  body  with 
respect  to  the  wall  sinusoid. 
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CHAPTER  ONE 


INTRODUCTION 

A  subnarine  running  parallel  to  a  wall  in  an  otherwise 
unbounded  fluid  experiences  forces  and  moments  which  are  a 
function  of  the  velocity  of  the  submarine  as  well  as  its 
proximity  to  the  wall.  These  forces  and  moments  arise  due  to 
the  reduced  pressure  between  the  hull  and  wall  caused  by  flow 
velocities  in  this  region  being  greater  than  anywhere  else 
around  the  body  surface;  this  phenomenon  is  commonly  referred 
to  as  the  Venturi  effect. 


The  force  developed  always  tends  to  pull  the  submarine 
towards  the  wall.  However,  the  direction  of  the  moment  is  a 
function  of  ocdy  shape.  Within  potential  theory,  no  net  moment 


act 

o  n 

a  body 

with  fore  and 

aft  symmetry 

near 

a  flat 

a  l  .  F  o  r 

a 

modern 

submarine 

with  blunt  bow 

and 

tapering 

. :  '  t  e  r  o  o  i  y  , 

the 

flow  over  the  bow 

i  s 

accelerated 

much 

more  than 

C  ; .  a  t  ever 

the 

3 1  e  r  n  , 

resul ting 

in 

a  moment  which 

tends  to 

rotate  the 

bow 

toward 

the  flat  wall. 

This 

problem  is 

of  practical 

interest  for 

a 

submarine 

operating  in  the  vicinity  of  the  ocean  bottom.  These  forces 
anc  moments  are  destabilizing  to  the  vessel’s  motion  and  will 
pull  the  subnarine  off  its  intended  course  into  the  bottom 
unless  accounted  for.  However,  if  the  force  and  moment  acting 
on  a  body  in  this  situation  can  be  predicted,  then  a  control 


system  can  be  devised  which  compensates  for  these  destabilizing 


t:.v_r>  :'or:e  calculations  approach  Hevman  *  s  theory  near  the 
wall,  nowever,  tney  rr.ove  closer  to  the  Lagally  results  as  the 
separation  distance  is  increased. 

All  tnree  techniques  predict  "bow-in"  yaw  moments  which  is 
physically  correct  for  a  blunt-nosed  modern  submarine  with 
tapering  afterbody.  Lagally  and  segmented  theory  calculated 
moments  are  comparable,  both  becoming  negligible  at  a 
separation  distance  of  approximately  three  to  four  hull 
diameters.  However,  results  from  both  of  these  methods 
differed  significantly  from  those  using  Newman's  theory. 

For  the  purpose  of  comparison,  the  normalized  drag 
(composed  primarily  of  frictional  drag  with  some  form  drag)  on 
this  unappended  modern  submarine  hull  form  in  a  real  but 
infinite  fluid  is  approximately  F 'drag  =  4.7.  It  is 
interesting  to  note  that  the  predicted  transverse  force  is  of 
the  same  order  of  magnitude. 

In  Figure  8,  force  calculations  are  compared  with  the 
results  of  model  tests  performed  on  an  actual  modern  submarine 
hull  form  in  the  vicinity  of  a  wall.  Although  the  tests  were 
done  usin^  an  appended  hull,  the  force  in  this  case  should  not 
differ  too  much  from  the  unappended  condition.  Lagally  and 
segmented  theory  force  calculations  are  shown  to  correlate  well 
with  the  test  results. 
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Figure  6.  The  Transverse  Force  Acting  on  a  Modern  Submarine 
Hull  Form  in  the  Vicinity  of  a  Flat  Wall. 


results 


As  a  function  of  separation  distance  S  (Figure  5),  the 
transverse  force  and  yaw  moment  on  a  representati ve  modern 
submarine  hull  form,  calculated  using  both  the  Lagally  and 
segmented  theories,  are  plotted  and  compared  with  Newman's 
slender  body  theory  in  Figures  6  and  7.  For  segmented  theory, 
a  two-dimensional  analysis  is  used  along  the  length  of  the  body 
except  for  the  region  extending  from  the  bow  aft  a  distance  of 
one-twent ieth  of  the  overall  length.  A  three-dimensional 
analysis  is  used  on  this  segment. 


FIGURE  5:  Separation  Distance  S 


Axial  force  was  calculated  using  Lagally's  theorem  to  be 
essentially  zero  as  anticipated  near  a  flat  wall. 

S  ranges  from  one-half  to  five  body  diameters. 

2  2  -4 

Force  is  normalized  by  1/2  pU  L  x  10  ;  moment  is  normalized  by 

2  3  -5 

1/2  pU  L  x  10  .  (U  is  body  forward  velocity;  L  is  body 

length.)  Primes  are  used  to  indicate  normalized  quantities, 
for  example  ,  F1  =  F/(l/2  pU2|_2  x  10*4). 

As  expected,  both  transverse  force  and  moment  increase 


rapidly  as  the  body  is  moved  closer  to  the  wall 


Segmented 


The  segmented  theory  presented  in  this  thesis  can  be  used 


to  calculate  only  the  transverse  force  on  a  body;  a  method  of 
determining  the  axial  force  is  still  to  be  developed.  However, 
since  no  axial  force  is  anticipated  for  the  flat  wall  case, 
this  poses  no  problem  at  this  time. 

Segmented  theory  offers  a  distinct  advantage  over 
Lagally's  theorem  in  that  the  appendages  on  a  body  can  be 
accounted  for  in  the  calculations.  By  merely  including  the 
appendages  when  calculating  the  added  mass  of  a  segment,  the 
force  and  moment  on  an  appended  hull  can  be  determined.  For 
the  purposes  of  this  thesis,  however,  the  appendages  will  not 
be  included  in  the  segmented  theory  calculations  in  order  to 
remain  consistent  with  the  other  methods  against  which  this 
technique  is  compared. 

Lagally's  theorem  predicts  the  moment  acting  on  a  body 
based  on  the  location  of  the  singularities  which  describe  its 
shape,  which  for  a  modern  submarine  hull  form  are  concentrated 
at  the  bow  and  stern.  In  actuality,  the  forces  which  generate 
the  yawing  moment  act  over  the  entire  body  length.  Because  the 
segmented  theory  approach  determines  the  transverse  force  at 
each  segment  along  the  body  length,  this  theory's  moment 
results  should  be  more  accurate  that  those  found  using 
Lagally's  theorem. 
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Figures  3  and  4  show  the  source  and  dipole  strengths 
during  the  iteration  process  of  sizing  a  singularity  system  to 
represent  a  modern  submarine  hull  fora.  (The  body's  velocity 
is  one  foot  per  second  and  is  located  with  a  separation 
distance  of  one-tenth  of  a  body  diameter  between  hull  and 
wall.)  With  iteration,  the  source  strength  remains  essentially 
constant,  reflecting  negligible  (as  compared  with  the  forward 
velocity)  induced  longitudinal  flow  along  the  body  and  image. 
Hov/ever,  dipole  strength  is  shown  to  increase  significantly 
between  the  first  and  fifth  iterations. 

As  the  body  is  moved  closer  to  its  image,  the  number  of 
iterations  until  convergence  increases  since  the  velocities 
induced  on  each  oody  are  greater.  For  this  reason,  an 
iteration  analysis  was  performed.  It  was  found  that  for 
separation  distances  down  to  one-half  of  a  diameter, 
convergence  is  guaranteed  within  three  iterations;  for 
distances  down  to  one-tenth  of  a  diameter,  five  iterations  are 
required;  for  distances  less  than  one-tenth  of  a  diameter,  the 
number  of  iterations  until  convergence  increases  rapidly.  In 
tue  "immediate"  vicinity  of  the  wall,  however,  the  iteration 
procedure  is  divergent  and  therefore,  the  velocity  potential  in 
this  situation  cannot  be  predicted  using  this  method. 

Once  the  velocity  potentials  for  the  body  and  its  image 
have  been  converged  upon,  the  force  and  moment  on  the  body  can 
be  determined  using  either  Lagally's  theorem  (Appendix  C)  or 
segmented  theory  (Appendix  D). 


A  body  in  the  proximity  of  a  flat  wall  can  therefore  be 
modeled  by  establishing  the  velocity  potential  for  two 
geometrically  similar  bodies  of  equal  size  moving  parallel  to 
one  another  at  constant  forward  velocity  in  an  otherwise 
infinite  fluid.  Appendix  B  presents  the  techniques  necessary 
to  establish  this  velocity  potential. 

CALCULATING  XH£  FORCE  MR  MOHBflT 

Computer  programs  (presented  in  Appendix  F)  to  calculate 
the  force  and  moment  on  a  modern  submarine  hull  near  a  wall 
were  created  based  on  the  following  considerations: 

The  body  and  its  image  can  both  be  sized  independently 
using  continuous  source  distributions  in  a  uniform  flow.  When 
these  two  distributions  are  brought  into  proximity,  each  will 
induce  velocities  over  the  other  body's  surface  which  will 
disturb  the  body  boundary  conditions  previously  satisfied. 
However,  by  using  a  slender  body  approximation,  a  continuous 
dipole  distribution  can  be  sized  to  restore  the  surface 
boundary  conditions. 

It  can  be  seen,  though,  that  these  boundary  conditions  are 
not  yet  fully  satisfied  since  the  velocities  induced  on  each 
body  are  now  caused  not  only  by  the  source  distribution,  but  by 
the  newly  introduced  doublet  system  as  well.  However,  the 
singularity  distributions  used  to  model  both  bodies  can  be 


resized  until  ultimately,  a  system  of  sources  and  dipoles  which 
satisfies  all  boundary  conditions  is  converged  upon  (Reference 


The  axisymmetric  slender  body  alone  in  an  Infinite  fluid 
can  be  modeled  using  a  line  distribution  of  sources  as 


C  J 


explained  in  Appendix  B,  "Modeling  an  Axisymmetric  Body." 

The  interaction  between  the  body  and  wall  is  taken  into 

account  by  meeting  the  boundary  condition  of  zero  flow  normal 

M.i  -  o 

to  the  wall,  that  is,  3n'c  .  This  is  done  using  the 

*wall 

method  of  images  in  which  the  wall  is  "replaced"  by  an  image 
(which  corresponds  to  the  body)  located  symmetrically  beyond 
the  wall.  Any  motion  of  the  body  is  "mirrored"  by  its  image 
and  therefore,  each  body  contributes  a  flow  normal  to  the  wall 
equal  in  strength  yet  opposite  in  direction.  The  resulting  net 


flow  across  the  wall  is  thus  zero. 


v*4>=o 
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FIGURE  2: 


Body  and  its  Image.  (Note  That  1^  =  0  Along  the  Line 

o  n 

of  Symmetry  Between  the  Two  Bodies) 


& 
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CHAPTER  TWO 


flat  MALL  amlxsxs 

The  force  and  noaent  acting  on  an  axisyianetric  body  moving 
parallel  to  a  flat  wall  can  be  determined  using  either 
Lagally's  theorem  (Appendix  C)  or  segmented  theory  (Appendix  D) 
after  a  velocity  potential  function  <P  has  been  constructed 
which  meets  the  appropriate  boundary  conditions.  The 
applicable  boundary  conditions  are: 

1.  Zero  flow  normal  to  the  body  surface, 


2.  Zero  flow  normal  to  the  wall, 


that  is,  u  _  3<p| 
vn  ”  3n 


Swall 


where  is  the  derivative  in  the  direction  of  the  unit  normal 


n  out  of  the  fluid.  The  governing  equation  throughout  the 

2 

fluid  domain  is  Laplace's  Equation,  V  <J>=0,  which  is  an 


expression  of  the  conservation  of  mass  for  a  potential 


function . 


FIGURE  1:  Boundary  Conditions  for  a  Body  Near  a  Flat  Wall 


the  wall 


Next,  the  irregular  wall  problem  will  be  addressed.  As  a 
first  step  towards  analyzing  the  case  of  an  irregular  wall,  a 
sinusoidal  wall,  varying  in  longitudinal  direction  only,  will 
be  investigated  to  simplify  the  problem.  The  problem  will  be 
limited  further  to  consider  the  "quasi-static"  case  in  which 
the  body  and  wall  are  fixed  in  a  uniform  flow  field  of  constant 
velocity  parallel  to  the  body  axis  at  infinity. 

The  force  and  moment  will  be  analyzed  over  separation 
distances  ranging  from  one- ha  If  to  five  body  diameters  between 
the  hull  and  wall.  At  distances  greater  than  five  diameters, 
it  is  anticipated  that  the  force  and  moment  acting  on  the  body 
will  be  negligible.  Distances  less  than  one-half  diameter  are 
not  considered  to  be  of  practical  importance  due  to 


beyond  the  wall 


Using  this  analysis,  an  approximate  velocity 


potential  was  ootained  from  which  the  pressure  distribution  on 
the  body  could  be  calculated.  The  force  and  moment  could  then 
be  determined  by  integrating  the  pressure  distribution  over  the 
body  surface.  However,  Eisenberg  failed  to  account  for 
velocities  induced  on  each  body  by  its  image  and  therefore, 
this  approximation  is  good  only  in  the  situation  where  the  body 

and  its  image  are  far  apart. 

2 

Newman  applied  slender  body  theory  to  a  set  of  images  and 
accounted  for  the  induced  velocities  by  offsetting  the  source 
distribution  an  appropriate  distance  from  the  body  axis.  This 
approach,  however,  is  only  good  for  bodies  in  the  immediate 
proximity  of  the  wall;  the  force  and  moment  fail  to  decrease 
as  rapidly  as  they  should  as  the  distance  between  the  wall  and 
body  increases. 


PRACTICAL  CONSIDERATIONS 

The  above-mentioned  approaches  are  concerned  with  bodies 
in  the  vicinity  of  a  flat  wall.  However,  the  ocean  floor  is 
generally  not  flat,  but  rather,  an  irregular  surface  which  can 
se  described  only  in  a  statistical  manner.  The  irregular  wall 
proolen,  therefore,  warrants  investigation. 

First,  the  case  of  an  axisymnetric  slender  body  moving 
with  constant  velocity  parallel  to  a  flat  wall  will  be 
approached  using  an  axial  distribution  of  singularities  to 


model  the  body  and  a  corresponding  image  system  to  account  for 


effects.  By  using  a  high  resolution  sonar  to  sense  the 
boundary  surface  contour,  such  a  control  system  might  permit 
high  speed  operation  in  the  vicinity  of  the  ocean  bottom. 

PROBLEM  FORMULATION 

In  order  to  determine  the  force  and  moment  acting  on  a 
submarine,  the  problem  must  first  be  idealized:  the  case  of  an 
unappended  slender  ( radius/ length  <<  1)  axisymmetric  body  in 
the  vicinity  of  a  wall  with  no  other  bodies  in  proximity  shall 
be  considered.  The  fluid  is  assumed  to  be  ideal  (inviscid)  and 
incompressible. 

The  particular  body  shape  analyzed  will  be  one 
representative  of  modern  submarine  hull  forms,  described  by  the 
following  characteristics : 

Length  Overall  /  Diameter  11 
Forebody  Length  /  Length  Overall  .17 
Forebody  Fullness  Factor  2 
Afterbody  Length  /  Length  Overall  .44 
After  Fullness  Factor  3 

Refer  to  Appendix  A,  "Hull  Geometry  Description  of  a 
Modern  Submarine"  for  further  detail. 

PREVIOUS  AND  RELATED  WORK 
1 

Eisenberg  considered  the  problem  of  a  spheroid  moving  in 
the  proximity  of  a  wall  using  a  source  distribution  to 


represent  the  body,  and  a  corresponding  image  source  system 


The  Transverse  Force  Acting  on  a  Modern  Submarine 
Hull  Form  in  the  Vicinity  of  a  Flat  Wall:  Model 
Test  Results  Compared  with  Theoretical  Calculations. 


CHAPTER  THREE 


WALL  ANALYSIS 


For  a  flat  wall,  the  boundary  condition  of  zero  normal 
flow  Is  satisfied  exactly  by  using  the  method  of  images. 
However,  this  technique  cannot  be  used  if  the  wall  is  not  flat. 

Consider  the  case  of  two  axisymmetric  slender  bodies 
alongside  one  another  as  shown  in  Figure  9  in  which  one  body  is 
very  "large"  in  comparison  with  the  other. 


Uniform 

Flow 


FIGURE  9:  Small  Body  in  Proximity  of  a  Large  Body 
Representing  a  Wall 


Boundary  conditions  over  each  body's  surface  can  be  met 
using  the  same  slender  body  approximations  discussed  in  the 
previous  chapter  to  size  source  and  dipole  distributions  along 
their  axes.  If  the  large  body  is  in  fact  "large  enough,"  then 
from  the  location  of  the  small  body,  it  will  appear  as  if  it  is 
a  wall.  As  the  diameter  of  the  large  body  is  increased,  the 
"wall"  will  flatten  out  as  shown  in  Figure  10.  In  the  limit  of 
an  infinite  diameter,  it  will  become  flat  in  cross  section. 
Thus,  by  varying  the  radius  of  the  large  body  sinusoidally  in 


By  restricting  the  analysis  to  sinusoids  with  wavelengths 


on  the  order  of  the  small  body's  length,  and  by  keeping  the 
ratio  of  the  sinusoidal  amplitude  to  wall  wavelength  "small," 
the  slender  body  assumptions  for  the  large  body  should  be 
satisfied.  Therefore,  an  axial  source  distribution  can  be 
sized  to  define  the  large  body's  shape  in  a  longitudinal  flow 
essentially  constant  along  its  length.  The  forward  and  after 
ends  of  the  sinusoidal  "wall"  body  are  considered  to  be  out  of 
proximity  of  the  small  body  so  they  have  no  influence  on  the 
flow  between  the  body  and  wall. 

A  doublet  distribution  is  superposed  on  the  source  system 
to  counter  the  transverse  velocity  induced  on  the  sinusoidal 
wall  by  the  body.  This  dipole  strength  is  proportional  to  the 
local  transverse  velocity. 

If  one  defines  a  "significant"  transverse  velocity  as  a 
local  velocity  which  is  greater  than,  say  one  percent  of  the 
maximum  induced  transverse  velocity  along  the  wall,  it  can  be 
seen  that  for  a  body  in  the  immediate  vicinity  of  the  wall, 
the  transverse  velocity  will  not  be  significant  along  a  length 
much  greater  than  the  body  itself  as  shown  in  Figure  12. 


Range  of  Significant 
Induced  Transverse  Flow 


FIGURE  12:  Transverse  Velocity  Induced  on  a  Wall 
by  a  Body  in  Close  Proximity 
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The  range  of  significant  transverse  flow  can  be  viewed  as 
the  "effective"  length  of  the  large  body  from  a  dipole  3izing 
standpoint  since  the  doublet  strength  is  essentially  zero 
outside  this  region.  This  effective  length  is  constant  for  a 
given  body  at  a  given  distance  fron  the  wall.  However,  as  the 
body  is  moved  away  from  the  wall,  the  effective  length 
increases  since  the  maximum  transverse  velocity  decreases  much 
more  rapidly  with  separation  distance  than  do  the  lesser 
transverse  velocities  some  longitudinal  distance  away. 

A  basic  assumption  in  sizing  the  dipole  distribution  is 
that  the  large  body  is  slender.  Therefore,  if  the  effective 
length  of  the  sinusoidal  wall  is  fixed,  then  its  diameter 
cannot  be  increased  without  bound.  Because  there  i3  no  clear 
limit  as  to  when  a  body  is  no  longer  slender,  a  scheme  had  to 
be  devised  to  allow  the  large  body’s  diameter  to  be  as  large  as 
possible  to  obtain  the  best  approximation  of  a  wall,  though 
still  remain  "slender." 

If  the  amplitude  of  the  sinusoidal  wall  is  zero,  then  it 
is,  in  fact,  a  large  cylinder  modeled  using  only  a  doublet 
distribution  in  a  transverse  flow  field.  In  performing 
calculations,  it  was  found  that  as  the  diameter  of  the  large 
cylinder  was  increased,  the  force  on  the  small  body  generally 
grew  until  it  reached  approximately  the  flat  wall  limit 
obtained  using  the  method  of  images.  Upon  increasing  this 
diameter  further,  the  force  actually  decreased,  indicating  that 
the  "effective"  body  was  no  longer  "slender."  By  using  the 
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largest  diameter  for  which  the  large  body  is  still  slender,  the 
best  approximation  of  a  wall  is  obtained.  Of  course,  this 
diameter  increases  as  the  body  is  moved  away  from  the  wall. 
Table  1  lists,  over  a  range  of  separation  distances  from  one- 
half  to  five  body  diameters,  the  approximate  maximum  "wall" 
diameter  for  which  the  "effective"  body  is  still  slender. 

Thus,  the  velocity  potential  for  the  case  of  a  body  in 
proximity  of  a  sinusoidal  wall  can  be  found  by  sizing 
singularity  distributions  for  two  axisymmetric  bodies.  Only 
the  quasistatic  case  is  to  be  analyzed,  that  is,  the  body  and 
"wall"  are  held  fixed  in  a  uniform  stream.  Therefore,  since 
both  are  axisymmetric  slender  bodies,  their  singularity  systems 
can  be  sized  using  the  same  method  discussed  in  the  previous 
chapter  and  presented  in  Appendix  B.  The  body  shape  is  that  of 
a  modern  submarine  hull  form;  the  "wall"  shape  is  that  of  a 
body  of  revolution  with  its  radius  varying  sinusoidally  about  a 
mean  radius  r0  in  the  longitudinal  direction  (Figures  11  and 
i 3 ) .  Appendix  E  describes  the  wall  shape  in  more  detail. 

In  the  case  of  a  sinusoidal  wall,  an  axial  force  on  the 
oody  is  expected,  unlike  the  flat  wall  situation.  This  axial 
force,  the  transverse  force,  and  yaw  moment  can  be  determined 
using  Lagally's  theorem  (Appendix  C)  once  the  velocity 
potential  has  been  established.  Similarly,  the  transverse 
force  and  yaw  moment  can  be  found  using  a  segmented  theory 
analysis  (Appendix  D). 


TABLE  1:  Maximum  Diameter  for  Which  The  Sinusoidal  "Wall" 
Body  is  Still  Slender 


SEPARATION  DISTANCE,  S  DIAMETER  OF  LARGE  "WALL"  BODY 


(in  terms  of  hull  diameter,  a) 


(in  terms  of  d) 
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In  Figures  14  and  15,  the  results  of  calculations 
predicting  the  transverse  force  and  consent  are  plotted  for  a 
modern  submarine  hull  form  near  a  flat  wall.  Results 
calculated  using  the  method  of  images  are  compared  against 
those  using  a  sinusoidal  wall  with  zero  amplitude  for  both  the 
Lagally  and  segmented  theory  techniques.  The  sinusoidal  wall 
results  are  in  agreement  with  those  calculated  using  the  method 
of  images,  the  sinusoidal  wall  predictions  generally  being 
slightly  less. 

Figures  16  through  21  present  the  results  calculated  using 
Lagally's  theorem  for  a  submarine  in  the  vicinity  of  a 
sinusoidal  wall.  Wall  sinusoidal  amplitudes,  Awall»  are  one“ 
tenth  and  three-tenths  of  a  body  diameter;  the  wavelength, 
Awa11*  in  all  cases  is  equal  to  body  length;  separation 
distances,  S,  are  spaced  with  one-half,  one,  and  three  body 
diameters  betv/een  the  hull  and  mean  wall  position.  Values  of 
force  and  moment  which  are  negligible  at  a  body  distance  of 
three  diameters  from  the  wall  generally  are  omitted  from  the 
,'raf  hs  presented. 

Forces  and  moments  are  calculated  with  body  position 
varied  longitudinally  while  parallel  to  the  sinusoidal  wall. 
The  longitudinal  position,  d  ^  ,  is  defined  as  the  axial  distance 
between  the  origin  of  the  submarine  body  (amidships)  and  the 
nodal  point  of  the  sinusoidal  wall  as  shown  in  Figure  13. 


A  Comparison  Between  the  Transverse  Force  Obtained  Using  the  Method  of  images 
and  that  Using  a  Sinusoidal  Wall  Analysis  with  Zero  Amplitude  for  a 
Modem  Submarine  Hull  Form  Traveling  Near  a  Rat  Wall 


Appendix  A 


Hull  Geometry  Description  of  A  Modern  Submarine 


A  modern,  unappended  submarine's  hull  shape  can  be 
described  most  simply  by  dividing  the  body  into  three  distinct 
geometries:  The  forebody  can  be  approximated  as  an  "ellipse" 
of  revolution,  the  parallel  midbody  as  a  cylinder,  and  the 
after  body  as  a  "parabola"  of  revolution.  Figure  A-1  indicates 
these  geometries. 


—  <t 

Afterbody  Parallel  Midbody  Forebody 

FIGURE  A-1.  Modern  Submarine  Hull  Geometry 


Mathematically,  the  radii  along  the  length  can  be 


described  as: 

Forebody:  r(xp) 


i/nf 

1 


(A.l  ) 


Parallel  midbody:  r(x) 
Afterbody:  r(xA) 


xfl  NA 
rn[l  -(—*-)  ] 


(A. 2) 
(A. 3) 


"N  "  and  "N  "  are  commonly  referred  to  as  fullness  factors 
F  A 

and  describe  the  degree  of  fullness  of  the  fore  and  after 


sections  respectively. 
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A  methodology  was  presented  for  predicting  the  forces  and 
moment  acting  on  a  submarine  near  a  sinusoidal  wall.  Results 
were  presented  only  for  the  case  of  a  sinusoidal  wall  with 
wavelength  equal  to  body  length.  Two  interesting  analyses  to 
be  investigated  further  are:  the  effect  of  decreasing  wall 
wavelength  and,  the  effect  of  varying  the  wall  radius  over  a 
wide  range  of  amplitudes.  However,  3ince  a  slender  body 
approximation  is  used,  it  would  be  necessary  to  first  determine 
how  far  the  limits  of  this  theory  extend  if  wavelength  is 
decreased  for  a  given  wall  amplitude  or  if  amplitude  is 
increased  for  a  given  wavelength. 

The  purpose  of  modeling  a  sinusoidal  wall  was  to  create  a 
basis  from  which  to  analyze  the  irregular  wall  problem. 
Perhaps  an  irregular  wall  can  be  approximated  under  certain 
circumstances  using  an  effective  equivalent  combination  of 
sinusoids.  However,  further  investigation  into  the  aspects  of 
the  irregular  wall  problem  are  required  to  determine  how  it 
might  be  modeled  best. 


Theorem  Results;  AwaU  =1/10  Body  Diameter,  5=  5/10  Body  Dlam 


Figure  25.  The  Moment  on  a  Modern  Submarine  Hull  Form  Traveling  Parallel 

to  the  Mean  Position  of  a  Sinusoidal  Wall  Using  Segmented  Theory 
AwqJI  =3/10  Body  Diameter 


3/10  Body  Diameter 


=1/10  Body  Diameter 


transverse  and  longitudinal  components.  The  axial  force,  shown 
in  Figures  20  and  21,  is  180°  out  of  phase  with  the  moment  and 
alternates  between  providing  forward  thrust  and  adding  drag 
force  to  the  body,  depending  on  the  longitudinal  location  of 
the  body  with  respect  to  the  sinusoidal  wall. 

The  magnitude  of  this  axial  force  is  of  the  same  order  as 
that  of  the  real  fluid  drag  (F^rag  =  •  7  )  indicating  that  at 

high  speed,  a  vessel  might  have  difficulty  maintaining  velocity 


constant  while  under  the  influence  of  such 
oscillating  axial  force. 


sizable 


Figures  22  -  25  show  the  transverse  force  and  yaw  moment 
which  result  from  segmented  theory  calculations.  Segmented 
theory  results  are  compared  with  Lagally's  theorem  results  in 


Figures  26  and  27  for  S  =  0.5d  and  A. 


0 . 1 d  .  Segmented 


theory  results  are  similar  to  those  calculated  using  Lagally's 
theorem  with  the  exception  that  the  transverse  force  mean  value 
along  the  wall  length  and  amplitude  are  approximately  3/2  times 
those  obtained  using  Lagally's  theorem.  Also,  the  mean  value 
and  amplitude  of  the  moment  are  approximately  twice  those 
determined  using  the  Lagally  analysis. 

As  previously  stated,  the  segmented  theory  calculations 
should  be  more  accurate  than  those  using  Lagally's  theorem. 
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Forces  and  moment  are  shown  to  vary  approximately 
sinusoidally  with  longitudinal  position  of  the  body  along  the 
wall.  As  the  amplitude  of  the  wall  is  increased  by  a  factor  of 
three,  from  one  to  three-tenths  of  a  body  diameter,  the 
amplitudes  of  the  transverse  attraction  force,  yaw  moment  and 
axial  force  increase  by  approximately  the  same  factor.  As 
expected,  forces  and  moment  decrease  rapidly  with  separation 
distance . 

Plotted  in  Figures  18  and  19  are  the  calculated  moments 
corresponding  to  the  forces  presented  in  Figures  16  and  17. 
The  moments  are  90°  out  of  phase  with  the  corresponding  forces. 
"Stern-in"  moments  are  developed  as  well  as  "bow-in"  moments. 

Hear  a  flat  wall,  a  "bow-in"  moment  is  always  predicted 
for  a  modern  submarine  hull  form  traveling  in  a  direction 
parallel  to  the  wall.  This  is  because  bow  and  stern  are 
equidistant  from  the  wall  surface  and  since  flow  is  accelerated 
much  more  at  the  bow  than  at  the  stern,  the  attraction  force 
forward  overshadows  the  force  aft,  resulting  in  a  "bow-in" 
moment.  However,  for  a  sinusoidal  wall,  if  the  wall  surface  is 


closer 

t  o 

the 

stern 

than  to  the 

bow,  the  attraction 

force 

acting 

on 

the 

after 

section  can 

become  greater  than 

that 

forward 

,  creating 

a  net 

"stern-in" 

moment . 

An  axial  force  is  expected  on  a  body  near  a  sinusoidal 
wall  since  the  wall  surface  is  not  parallel  to  the  direction  of 
body  motion.  The  attraction  force  acting  on  the  body  due  to  a 


differential  length  of  wall  surface  can  be  decomposed  into 


o  Modam  Submortto  Hull  Form  Trawling  PoralM 
mi  of  a  SkwooMai  MMf  Uatig  Logall/o  Thooram 


The  slopes  of  the  hull  along  the  length  of  each 
can  be  determined  by  differentiating  equations  (A.1)  to 
with  respect  to  the  appropriate  X  coordinate: 
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Appendix  B 


Modeling  An  Axi symmetric  Slender  Body 

The  flow  around  a  slender  body  of  revolution  moving  with 
forward  velocity  II  in  an  infinite  ideal  fluid  can  be 
approximated  using  a  continuous  axial  source  distribution 
superimposed  on  the  velocity  of  a  uniform  stream. 


o)  Body  Travaftig  at  Voleclty  U 

FIGURE  B- 1  . 


b)  Find  Body  to  a  Uniform  Stroan 


If  the  body  geometry  is  described  as: 

r  =  f (x) 

where  r  is  the  body  radius,  then  the  combined  velocity 
potential  of  source  distribution  and  free  stream  is: 


<j>  ( x  )  =  -Ux  + 


jeIlL 


C  (X  -er  +  r‘] 


T/2 


d£ 


(  B  .  1  ) 


where  U  is  the  free  stream  velocity,  a  and  b  are  the  body 
endpoints  as  shown  in  Figure  B-1 ,  and  m( £)  is  the  source 
strength  at  an  axial  position  represented  by  the  "dummy" 
variable  £.  The  body  oriented  orthogonal  axes  x,y  and  z  are 


defined  in  Figure  B-1  with  the  origin  at  the  midships  position 


along  the  centerline  of  the  body 


The  radial  velocity  with  respect  to  the  body  oriented 


q r a d i a  1  Gan  *3e  determined  to  be: 


q  (x )  =  -  =  f* 

qradi al  '  '  dr  j . 


m(g)r  dg 
C(x  -  5)2  +  r2] 


which  can  be  rewritten  as: 


^radi al  ^  x  ^  r 


b-x ) 


'(a-x  ) 


[(^)2  +  1] 


(B.2) 


372  d(^~)  ( B  . 3  ) 


For  slender  bodies  of  which  <<  1,  where  L  is  the  body  length, 
the  term  [-£—  ~-X]2  >>  1  except  in  the  vicinity  of  £  =  x. 
Therefore,  m(  5)  can  be  replaced  by  ni(x)  and  the  integral  limits 


extended  to  -  00  to  +00  .  Substituting  n  s 
can  be  rewritten  as: 


~  x  ,  equation  (3.3) 
r 


"radial*10  = 


..  m(x) 


an  _  mUi 
-  (n2  +  1)3/2  "  r 


(  B  .  4  ) 


Since  the  body  is  slender,  a  ,  (x)  can  also  be 

Mradi al 


approximated  to  first  order  as 


d  r  a  d  i  a  1  ^  x  ^  ^ 


(  B  .  5  ) 


because  the  velocity  in  the  x  direction  is  essentially  U  along 
the  body  side. 

Equating  equation  ( B . 4 )  with  (D.5), 


Therefore,  the  source  strength  must  be: 


m(x)  =  Ur  ( x )  (B.7) 

An  axisynmetric  body  moving  steadily  along  its  axis  in  an 
infinite  inviscid  and  irrotational  fluid  experiences  absolutely 
no  hydrodynamic  forces.  However,  if  singularities  external  to 
the  bocy  are  introduced  into  the  fluid  field,  then  forces  will 
act  on  the  body.  These  forces  can  be  determined  using 
Lagally's  theorem  or  segmented  theory  if  the  internal  and 
external  singularities  are  known. 

When  the  axial  distribution  of  sources  is  brought  into  the 
proximity  of  other  singularities  external  to  the  body  of 
revolution,  velocities  will  be  induced  on  the  body  surface 
which  will  disturb  the  satisfaction  of  the  surface  boundary 
conditions  previously  employed.  In  order  to  re-establish  the 
body  surface  boundary  conditions,  it  is  usually  necessary  to 
introduce  additional  singularity  distributions  into  the  body 
which  are  images  of  the  externally  induced  flow.  These  image 
distrioutions  can  be  sized  exactly  only  for  spheres.  However, 
for  slender  bodies,  an  approximate  singularity  distribution  can 
oe  determined. 

Consider  the  case  of  two  geometrically  similar  bodies  of 
revolution  of  equal  size  moving  parallel  to  one  another  along 


FIGURE  B-2  Two  Geometrically  Similar  Bodies  of  Revolution 
Moving  Parallel  to  One  Another. 

Each  body  alone  in  an  infinite  fluid  can  be  modeled  as  a 
continuous  source  distribution  using  equation  (B.7),  however, 
when  two  bodies  are  placed  in  proximity,  velocities  induced  on 
each  body  by  the  other  will  distort  their  shapes. 

Consider  now  the  transverse  velocity  induced  on  one  body. 
Since  this  body  is  slender,  it  is  reasonable  to  assume  that 
over  any  short  axial  length,  the  transverse  velocity  which 
would  exist  if  the  body  were  not  present  is  approximately 
uniform.  An  axial  distribution  of  doublets  can  be  sized  which 
will  counter  this  cross  flow  velocity  and  restore  satisfaction 
of  the  boundary  conditions. 

Since  this  analysis  is  restricted  to  slender  bodies,  it 
can  be  further  assumed  that  the  body  will  be  cylindrical  over 
any  small  change  in  axial  location.  Therefore,  the  image  axial 
dipole  distribution  oan  be  sized  using  two  dimensional  theory. 

For  a  cylinder  of  infinite  length  in  a  uniform  transverse 
flow,  its  velocity  potential  may  be  expressed  as: 


*<x)  ■■  L 


(B  .8 ) 


where  5  is  the  strength  of  the  axial  distribution  of  doublets, 

constant  along  the  infinite  cylinder's  length;  y  and  r  are, 
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respectively,  the  transverse  and  radial  distances  from  the 
axis.  This  expression  is  valid  if 


?  = 


(  B  .  9  ) 


where  R  is  the  radius  of  the  cylinder  and  V  is  the  uniform, 
transverse  velocity. 


Treating  the  slender  axisymmetric  body  as  a  piecewise 
summation  of  cylindrical  components,  the  image  doublet 
distribution  strength  will  be: 

C(x)  =  r2(x)v(x)  ( B  .  1  0  ) 

where  jr(x),r(x)  and  V(x)  now  vary  along  the  body  length.  This 
approximation  is  good  provided  the  slope  of  the  body  is  not  too 
large;  near  the  blunt  forward  end  of  a  modern  submarine,  this 
approximation  is  very  poor. 

The  situation  of  two  bodies  of  equal  size  moving  parallel 
to  one  another  as  shown  in  Figure  B-2  can  be  considered  as  the 
case  of  a  body  and  its  image.  In  this  situation,  due  to  the 
symmetry  of  flow  about  the  x-y  plane,  the  dipole  distribution 
oriented  in  the  z  direction,  (x),  will  be  zero  along  the 
entire  body  length. 

The  induced  longitudinal  velocities  u(x)  can  be  accounted 
for  by  resizing  the  axial  source  strengths  accordingly,  that 
is, 

■.(»)  .  -  lu(?>  ■  U>  r(  x)  -j}lM  ( B .  1 1  ) 

However,  induced  longitudinal  velocities  will  generally  be 
insignificant  when  compared  with  the  free  stream  velocity. 
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Shifting  this  analysis  to  the  image  body,  it  can  be  seen 
that  the  image  has  velocities  induced  on  it  which  are  caused  by 
the  resized  source  and  newly  introduced  doublet  distributions 
of  the  body  which  will,  in  turn,  disturb  the  image  body 
boundary  conditions.  However,  by  iterating  this  procedure  of 
resizing  the  singularities  on  each  body,  convergence  can  be 
obtained  and,  ultimately,  all  boundary  conditions  satisfied 
(Reference  3 ) • 


Appendix  C 

Determining  the  Force  and  Moment  on  An 
Axisyrametric  Body  Using  Lagally's  Theorem 


The  Lagally  Force  and  Moment 


After  having  obtained  the  velocity  potential  description 
for  a  body  in  th^  vicinity  of  a  wall,  the  pressure  distribution 
around  the  body  can  be  obtained  using  Bernoulli's  Theorem.  By 
integrating  this  pressure  over  the  body  surface,  the  force  and 
moment  acting  on  the  body  can  also  be  determined. 

However,  if  the  body  is  modeled  as  an  axial  distribution 
of  discrete  sources  and  doublets,  the  force  and  moment  acting 
on  each  singularity  will  be  a  function  of  the  singularity  type 
and  strength  as  well  as  the  velocity  induced  at  its  location  by 
all  singularities  external  to  the  body.  By  summing  up  the 
force  and  moment  due  to  each  singularity  within  the  body,  the 
total  force  and  moment  can  be  determined.  Since  in  fact  the 
body  is  represented  by  a  continuous  axial  distribution  of 
sources  and  doublets,  the  summation  becomes  an  integration  of 
differential  forces  and  moments  along  the  axis. 

Lagally  was  the  first  to  reason  this  approach  and  derived 
simple  expressions  which  relate  the  force  and  moment  on  a  body 
to  the  singularity  distribution  which  describes  it.  The 
Lagally  force  and  moment  are  limited  to  steady-flow  problem 
analyses,  however,  this  is  consistent  with  the  constant 
velocity  situation  studied  in  this  thesis. 
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The  steady-state  Legally  forces  and  momenta  ban  be  expressed 

as: 

T(C)  -  -  4irp[m(OqU)  +  (cU)  •  V)q(U]  (C.l) 

HU)  *  7(0  x  fU)  +  47rp(qU)  x  ?(£))  (C .2 ) 

vhere  T(£)  and  M(C)  ****  the  differential  force  and  oooent, 

respectively,  acting  on  a  source  of  strength  m(  £ )  and  doublet 
of  strength'll;)  at  the  axial  location  described  by  the  position 

A  A  A 

vector  rt  5  )  *  5  1  ♦  0  j  ♦  ok  .  q(£)  i3  the 

velocity  induced  at  the  location  T*(£  )  by  all  singularities 
external  to  the  body  and  p  is  the  vater  density. 

Decomposing  equations  (C.l)  and  (C.2)  into  orthogonal 
components: 

Fx({>  *  -  4itp[m({)  q^U)  ♦  Uy(E)  fy)  qx Ct)3 
Fy(E)  -  -  4*p[*U)  qy(E)  +  ( ty  ( e )  'ly)  qy(t>3 
FjU)  *  -  4.p[m(C)  qr(e)  ♦  Uy(E)  -|y)qz(e)3 

(C  .3 ) 

MXU)  c-4tip£  cy(C)  qz(5)] 

MyU)  =  -  Ztz(K)  (C  .4 ) 

MZ(C)  =  CFy  (  C  )  +  47ipcy  ( C)  qx(U 

where  the  subscripts  x,  y,  and  z  denote  the  respective 
component  of  force  or  moment.  Analyzing  for  the  case  of  two 
axisymmetric  bodies  moving  parallel  to  one  another  in  an 
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otherwise  Infinite  fluid,  the  doublet  distributions  In  the  x  and 
z  directions,  and  Cz  are  both  zero  along  the  centerlines  of 
both  bodies.  Therefore,  these  terms  have  been  omitted  from 
equations  C.3  and  C.4. 

By  symmetry  of  the  induced  flow  about  the  plane  z  =  0, 

qz(  £ )  is  also  zero.  Therefore,  equations  (C.3)  and 

(C.4)  can  be  reduced  to: 

FXU)  -  -  4itp[m(C)  qx(C)  ♦  tyU)  fy  qxU)3 

ryU)  *  -  4ttpOU)  qy(C)  ♦  Cy(U  |y  qy(c)3  (C.3a) 

f2(c)  -  0 


MX(C)  -  0 

My(U  -  o  (C.4a) 

MZ(C)  *  ayU)  ♦  4irpcy(C)  qx(4) 


Induced  Velocities  and  Velocity  Gradients 


Next,  in  order  to  determine  the  force  and  moment  on  the 
body,  it  is  necessary  to  develop  expressions  for  the  velocities 
and  velocity  gradients  induced  on  the  body  by  external 
singularities . 

Consider  the  case  at  hand  in  which  two  bodies  of 
revolution  are  moving  along  their  axes,  each  parallel  to  the 
other,  at  constant  forward  velocity  U,  as  shown  in  Figure  C-1. 


>yyy;  v-v 


V-W. --V-V^Y -v 


FIGURE  C-1.  Two  Bodies  of  Revolution  Moving  Parallel  To  One 

Another  at  Constant  Velo.city  U. 

If  the  singularity  distribution  in  body  B1  is  composed  of  an 
axial  source  distribution  superimposed  on  a  transverse  dipole 
distribution  oriented  In  the  y  direction,  the  potential 
function  at  a  point  (x^,  y2  *  *2*  located  ln  B2  due  to  the 
singularities  at  point  (x^,  y^ ,  )  in  B1  can  be  described  as: 

—  m ( x - ,y «  ,Zj  ) 

♦(*2»yp»22^  = 

rl  2 

where  m(x„,  y, ,  z, )  and  r  (x„,  y, ,  z„)  are  the  source  strength 
111  sy  1  1  1 

and  dipole  strength  respectively  at  the  point  (x^,  y^  ,  z^>  and 
r12  *  /(x2  -  +  (y  ^  -  y-j  )2  +  (z2  -  z^)'d . 

For  the  problem  at  hand,  it  is  convenient  and  appropriate  to 
take  y^  =0  and  z^  =  0  since  the  singularity  distribution  is 
along  the  axis,  and  substitute  the  "dummy"  variable  £  for  x^  . 
In  order  to  find  the  total  velocity  potential  at  a  point 
( x  2 »  y 2*  ^  )  due  to  only  the  axial  singularity  distributions  in 


cy(xT .yj )  (y2  -  ) 


( C  .  5 ) 
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body  B1 ,  it  is  necessary  to  integrate  equation  (C.5)  over  the 
length  of  B1 . 


That  is, 


fb  t 

^  (  x2  ’^2  ’  Z2  ^  ~  |" 


?y(c)  y2 


(  C  .  6  ) 


where  r]2  =  [(x2  -  i)2  +  y  2  +  z2} 

This  integration  can  be  performed  numerically  using  Simpson's 


Rule. 


Velocities  induced  at  points  located  in  B2  can  be 
determined  by  differentiating  the  velocity  potential  and 
evaluating  it  at  the  desired  location. 

„  .  3«*.(xz ,y2.z2)  fb  j<V5)«U)  _  3(x2-£);y(c)y2| 

b  (C.7) 

-  i.  ..  .  ,  .  11 f  .  3VV5))  .. 


qy(x2’y2,Z2)  = 


( C  .  8 ) 


Similarly,  the  velocity  gradients  ~( q  )  and  -r— (q  )  can 

9y  x  dy  y 

be  obtained  by  differentiating  the  velocities: 

3  ,,  *  (  3{x,-£)[c v(5)-*(£)y 2  .]  1 5 ( x.-E )y/  c  ( 5 )  ) 

=  j— - h - - - 2 - 7  -  (dS 

a  v  r  1 2  r12  7 

( C  .  9  ) 

..  ,  m  -  f  /ni( C }  .  3y2[3S(C)"y 15 y?3Z  (t)  ) 


-^Cq  (x2  ,y2,22)  3  ■ 


a  ir . 


'  12  1 2  1  12  ' 

(C.  10) 

Expressions  for  the  differential  force  and  moment  acting 
along  the  body  length  can  now  be  determined  by  evaluating 
equations  (C.7-10)  at  a  point  (  £2>0,0  )  along  the  B2  axis,  and 
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substituting  the  results  into  equations  (C.3a)  and  (C.Ua). 
(MOTE:  C-|  and  ?2  are  introduced  here  to  distinguish  between 
the  dummy  variable  £  on  body  B1  and  on  B2.)  The  total  force 
and  moment  acting  on  B2  can  be  determined  by  integrating 
equations  (C.3a)  and  (C.4a)  along  its  length. 


Appendix  D 


The  Force  and  Moment  on  an  Axisymnetric  Body 
Using  Segmented  Theory 


If  a  slender  axisymmetric  body  is  maintained  at  constant 
velocity  in  an  accelerating  fluid,  the  lateral  force  on  the 
body  can  be  approximated  using  a  "segmented"  theory  proposed  by 
Abkowitz.  "Segmented"  refers  to  dividing  the  body  into 
vertical  segments  and  calculating  the  force  on  each  segment 
using  either  a  two  or  three-dimensional  flow  analysis, 
whichever  is  appropriate.  In  the  two-dimensional  case,  the 
segments  correspond  to  "strips"  used  in  strip  theory. 

Along  most  of  the  body's  length,  it  can  be  assumed  that 
the  flow  is  locally  two-dimensional.  Therefore,  by  dividing 
this  portion  of  the  body  into  a  series  of  thin  transverse 
slices,  the  force  on  each  segment  can  be  found  and  integrated 
along  the  body's  length  to  determine  the  total  force.  This 
method  of  analysis  assumes  that  the  flow  at  any  segment  is 
independent  of  the  flow  at  any  other  location,  that  is,  there 
are  no  hydrodynamic  interactions  with  any  other  segments  along 
the  body. 

However,  near  the  blunt  bow  of  a  modern  submarine,  three 
dimensional  end  effects  are  significant  and  must  be  taken  into 
consideration.  These  end  effects  can  be  accounted  for  by  using 
a  three-dimensional  analysis  which  treats  the  bow  segment  as  a 
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KUMSTA  -  Length  of  the  ellipsoidal  bow  (in  number  of  stations 
from  the  forward  perpendicular). 

PHIXY  -  ^xy 

mil  -  ^yy 

H  -  Distance  between  a  singularity  point  and  a  field  point. 

RAD  -  Local  radius  of  the  body. 

KADP  -  Local  slope  of  the  body. 

RO  -  Fluid  density. 

£  -  Station  spacing  of  the  body. 

SZ  -  Station  spacing  of  the  ellipsoidal  bow. 

SOURCE  -  Local  source  strength  along  the  body  axis. 

SOURCE1./ AL  -  Local  source  strength  of  the  "large"  body  used  to 
approximate  a  sinusoidal  wall. 

SU  -  Station  spacing  of  the  sinusoidal  wall. 

THETA  -  Phase  of  the  sinusoidal  wall  as  seen  from  the  midships 
position  of  the  body  in  proximity. 

JI  -  Induced  longitudinal  flow. 

X  -  Induced  transverse  flow. 

VI IiTHASSFOR  -  Sum  of  the  displaced  mass  and  added  mass,  that 

is,  the  virtual  mass  of  the  ellipsoidal  bow. 

V 0 L F 0 R  -  Displaced  volume  of  the  ellipsoidal  bow. 

VJALRAD  -  Local  radius  of  the  large  body  used  to  approximate  the 
sinusoidal  wall. 

'.I ALRADP  -  Local  slope  of  the  large  body  used  to  approximate  the 
sinusoidal  wall. 

X  -  Distance  between  the  centerlines  of  two  bodies  in 


proximity. 


Listing  of  Significant  Variables  used  in  The  Computer  Programs 


ADMASSF ACT  -  Ratio  of  the  added  mass  to  the  displaced  mass  of 

the  ellipsoidal  bow;  this  factor  can  be 

4 

determined  using  Article  115  of  Lamb  . 

AMP  -  Amplitude  of  the  sinusoidal  wall. 

AVG ACC  -  Average  fluid  acceleration  acting  on  the  ellipsoidal 


3D  -  Ratio  of  large  sinusoidal  body  diameter  to  the  diameter  of 
the  suomarine  hull  form  in  proximity. 

CENTVOL  -  Longitudinal  center  of  volume  of  the  ellipsoidal  bow 
(from  the  forward  perpendicular). 

PI  A  -  Maximum  diameter  of  the  submarine  hull  form. 

DIPOLE  -  Local  dipole  strength  along  the  body  axis. 

DIPQLEVf AL  -  Local  dipole  strength  of  the  large  body  used  to 
approximate  a  wall. 

FVEL  -  Forward  velocity  of  the  body  in  ft/sec. 

FXTOTAL  -  Total  force  acting  on  the  body  in  the  X  direction. 
FYTOT AL  -  Total  attraction  force  on  the  body. 

2tt 

KWALL  -  Wave  number  of  the  sinusoidal  wall;  XWALL  =  , - 

Awall 

LAE  -  After  body  length  of  the  modern  submarine  hull  form. 

LENGTH  -  Length  (overall)  of  the  body. 


LFB  -  Forebody  length  of  the  modern  submarine  hull  form. 
Il.il’:.  I'M  -  Simpson's  multipliers. 

KASSFOR  -  Fluid  mass  displaced  by  the  ellipsoidal  bov/. 
iiZTQTAL  -  Total  yaw  moment  acting  on  the  body. 


-  Number  of  stations  along  the  sinusoidal  wall 


Appendix  E 

Sinusoidal  Wall  Geometry 


Figure  E.1  defines  the  sinusoidal  wall  geometry  used  in 
this  thesis.  The  wall  is  modeled  using  a  large  axisymmetric 
body  whose  radius  varies  sinusoidally  in  the  longitudinal 
direction  about  a  mean  radius  r^ .  Mathematically,  the  wall 
surface  can  be  described  as 


n<xw>  =  r0  "Awa11  SINtKwa11(xw  -  dL  >] 


where  A  ,  ,  is  the  sinusoidal  amplitude;  K  ..is  the  wave 
wall  wall 

number  of  the  wall,  that  is  K  .  „  =2ir/A  where  1  is  the  wall 

wall  wall  wall 

wavelength.  4^  is  the  longitudinal  distance  from  the  origin  of 
the  body  centered  orthogonal  axes  (x=0)  to  the  origin  of  the 
wall  sinusoid  (x^sO)  and  expresses  the  phase  of  the  sinusoid  at 
a  point  on  the  wall  located  transversely  across  from  the  body 


origin . 


— 1 - 1 


FIGURE  E.1. 


Sinusoidal  Wall  Geometry 


Therefore , 


[vtO]  - 


p  u  +  V 
3x  3y 


3^4)  3  4>  3  ^ (ft  3  (p 

3x3 y  3x  3y2  gy 


Because 


(D.6) 

(D.7) 


,  and  since  M  >>  -JL£ 

3y 

and  3  <|> 

3  xgy 


is  of  the  same  order  as  3  9  , 

ay2 

the  contribution  of  the  y  component  of  transverse  acceleration 
is  negligible  compared  with  the  longitudinal  component.  For 
completeness,  however,  all  terms  will  be  retained. 


ybMi* 

^  vWV  V  '/  \ 

-  -  >'• *  -  *  «  '  .  - 
-  V  - A  >  "•  >  ^  ’ 

k“*  ’.  v  •  *’•  W  *  V**  * 

dividing  by  two  the  value  calculated  to  be  the  added  mass  of 

the  equivalent  ellipsoid  of  which  the  he mi-ellipsoid  is  half. 

4 

Lamb  ,  in  article  15,  presents  the  information  necessary  to 
determine  the  added  mass  of  an  ellipsoid. 

The  moment  on  the  forebody  will  be  the  product  of  this 
force  and  its  moment  arm  from  the  origin.  The  center  of  force 
of  the  hemi-ellipsoid  is  approximated  to  be  at  its  center  of 
volume . 

For  the  analyses  herein,  the  Froude-Krylov  assumption 
shall  be  invoked,  that  is,  it  will  be  assumed  that  the  body‘s 
presence  does  not  affect  the  pressure  or  acceleration  fields 
around  it.  The  acceleration  can  therefore  be  determined  from 
the  local  velocity  potential  generated  by  all  singularities 
external  to  the  body  as  well  as  any  uniform  flow  in  which  the 
body  is  located. 

The  substantial  acceleration  can  be  expressed  as: 


9  v  ,  3  v  3x  3v  3y  3v  3z 

3 1  3x  3t  3y  9t  3z  3t  • 


(D.5) 


Because  the  problem  analyzed  is  for  steady  flow,  =  o  • 

a  t 

Also,  since  the  acceleration  is  to  be  calculated  along  the  body 


axis,  by  symmetry  of  the  flow  induced  by  the  other  body, 


-70- 


force  on  a  se^nent  due  to  its  presence  in  an  acceleration  field 
which  causes  a  pressure  gradient  at  the  segment's  location. 
This  force  is  dependent  on  the  volume  of  the  segment  and  is 
generally  referred  to  as  the  Froud e-Krylov  force.  The 
^22  Ft  Cv(£)3  terra  represents  the  effect  of  fluid 

acceleration  relative  to  the  segment.  Thi3  is  the  so-called 
added  mass  force  which  depends  on  the  segment  shape. 

Since  the  added  mass  of  a  two-dimensional  cylinder  is 

p 

equal  to  pA  and  A  is  equal  to  Tir(C)  *  where  r(£)  is 
the  local  body  radius,  equation  (D.1)  can  be  rewritten  as: 

FUteral<«>  =  ^U)2  St  <°'2> 

The  total  hydrodynamic  force  on  the  after  segment  can  now  be 
determined  by  integrating  equation  (D.2)  along  the  segment'3 
length : 

rc 

lateral  '  2lt»  r<5>2  ST[v<5)]  de  (D'3) 

J  31 

Similarly,  the  yaw  moment  on  the  after  section  can  be  found: 
fC 

Mz  =  2frp  5rU)2  {^OU)]  d?  (0.4) 

d 

The  force  on  the  heni-ellipsoid  can  be  obtained  by  taking 
the  product  of  the  average  local  substantial  acceleration  and 
the  sum  of  its  three-dimensional  added  mass  plus  the  displaced 
fluid  mass.  This  force  represents  the  sum  of  the  Froude-Krylov 
and  added  mass  forces  as  previously  discussed  with  the 
exception  that  a  three-dimensional  analysis  applies  here.  The 
added  mass  of  the  hem.i-ellipsoid  can  be  approximated  by 


x  =  c  X  =  b 


FIGURE  D-1.  Slender  Axisymmetric  Body  Divided  Into  Two  Sections 
An  Afterbody  to  be  Analyzed  Using  Two-Dimensional 
Theory  and  a  Hemi-Elllpsoidal  Bow  to  be  Analyzed 
Using  Three-Dimensional  Theory. 

The  force  acting  on  a  two-dimensional  segment  can  be 
expressed  as: 

FUteral<5>  ■  <  OA  *  A22)  {fcCvU)]  (D.l) 

where  £  is  a  dummy  variable  representing  the  axial  location  of 
the  segment;  p  is  the  fluid  density;  A  is  the  cross  sectional 
area  of  the  body;  (  p  A  is  the  local  mass  of  fluid  displaced  per 
unit  length);  A  ^  is  the  local  transverse  added  mass  of  a  two- 
dimensional  cylinder  with  cross  sectional  area  A.  [v(£)] 

is  the  substantial  derivative  of  the  local  transverse  fluid 
velocity,  that  is,  an  expression  of  the  transverse  fluid 
acceleration  as  viewed  from  a  coordinate  system  moving  with  a 
fluid  particle. 
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